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Abstract 

We present a realistic model for transferring the squeezing or the entanglement of optical field 
modes to the collective ground state nuclear spin of 3 He using metastability exchange col·lisions. 
We discuss in detail the requirements for obtaining good quantum state transfer emciency and 
study the possibility to readout the nuclear spin state optically. 
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I. INTRODUCTION 



Helium 3 àtoms in their ground state possess a purely nuclear spin / = 1/2. Such 
spins are well-isolated from the environment and show extremely long coherence times. 
jongitudinal coherence times 7\ of several days are measured in room temperature samples 
ÏJ. The transverse coherence time T2, which would be as long as T\ in zero magnètic 
field, is usually limited by magnètic field inhomogeneity if no special precaution is taken. 
Transverse coherence times of several hours are observed in a very low field |2j|. These 
very long coherence times originate from the weakness of magnètic coupling on the one 
hand, and from the absence of electrical coupling on the other hand, as there is no elèctric 
quadrupole coupling within the ground state for spins 1/2. It is tempting to exploit such 
long-lived coherence for quantum information purposes. In a previous letter we studied 
the possibility to transfer the squeezing of a cavity mode to 3 He nuclear spins. We showed 
that the squeezeing could be stored and retrieved from the àtoms, thus realizing a quantum 
memory □□□□□□ . For the sake of simplicity we presented in our letter a simplified 
model involving only two sublevels in the metastable state and gave numerical results for 
the more complicated case of 3 He. In this paper we concentrate on 3 He and treat in detail 
this more realistic case. 

Section [iï] is devoted to metastability exchange col·lisions. We derive linearized 
Heisenberg-Langevin equations describing the exchange collisions from the standpoint of 
quantum fluctuations. In section IIHI we describe the model for squeezing transfer from a 
squeezed vacuum mode of the electromagnètic field to the àtoms. Numerical results are 
shown and discussed in section ITV1 In section 0we obtain analytical results in the adiabatic 
elimination limit for the optical coherences and the cavity field. Section IVll is devoted to the 
readout scheme of the quantum memory. In section IVIU as a straightforward application 
of our scheme, we consider the possibility of creating long-lived quantum correlations be- 
tween two macroscopic spins, in the move of the successful experiment in Copenhagen |ïflj| . 
in which two macroscopic spins were entangled for 0.5 ms, but on a completely different 
timescale. Finally, in section IV1111 we use a toy model to explore the consequences of an 
imperfect polarization of the àtoms on our squeezing transfer scheme. 
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II. METASTABILITY EXCHANGE COL·LISIONS IN HELIUM 3 



Over forty years ago, Colegrove, Schearer and Walters demonstrated a technique 
to polarize 3 He relying on (i) an optical interaction on an infrared transition from the 
metastable 2 3 S triplet state to the 2 3 P triplet state, and (ii) metastability exchange collisions 
between àtoms in the ground state and in the metastable state. During such a collision, two 
àtoms exchange their electrònic degrees of freedom so that the metastable atom, oriented by 
optical pumping and with a nuclear polarization due to hyperfme coupling in the metastable 



state, becomes a polarized ground state atom 12]. This technique called metastability 
exchange optical pumping is currently used to prepare polarized samples for nuclear physics 
experiments as well as in nuclear magnètic resonance imaging applications [ï^ . 

In what follows we suggest that metastability exchange collisions can also be used to 
transfer quantum correlations to the ground state nuclear spin of 3 He. 



A. Equations for the one-body density matrix elements 

Partridge and Series describe a metastability exchange (ME) collision in terms of the 
one-body density matrices representing the infernal states of two colliding àtoms that we 
name pjf and p^ for the ground and metastable state, respectively. The density matrices 
after the collision, and , are given by 

Tr rt at 

e Pm (1) 
pf <g> Tr n p^ 

where Tr e and Tr n are trace operators over the electrònic and nuclear variables. 

Let us consider n metastable and N ground state independent àtoms. We introduce 



at' 



P'l 

Pm 



Npf, p m 



np^, and the same for p' and p' m . 




To represent the state of the system, we will use the density matrix p defined by: 

(2) 

V P%) 

Note that Trp = n + N and that we neglect all coherences between the ground and the 
metastable states. The matrices p g and p m evolve according to 

-7/Pg + 7/Pg 

TmPm Tm P m 



dt^S 



(3) 



where 7/ and 7 m are the metastability exchange collision rates in the ground and metastable 
states respectively 

7m ^^fexc 7/ ^7exc (4) 

with 7 exc a rate depending on the metastability exchange cross section, the relative velocity 
of the àtoms and the volume explored by the àtoms. 

The calculation of dp/dt is performed by expressing p in the decoupled spin basis of the 
nuclear spin I = \ and the total electrònic spin J = S = 1 in the metastable state, followed 
by a projection onto the hyperfme states (eigenstates of the total momentum operator F 
and F z ) which we name from 1 to 6 as in figured The explicit evolution equations for the 
density matrix elements are given in the Appendix. The fully polarized state in which all the 
àtoms are in the sublevel with highest angular momentum projection along z is stationary 
for equations (jïïj). 




m F =-| mp= 2 mp= 2 mp= 2 

FIG. 1: Sublevels 1 to 6 are metastable; 9 and are the ground state sublevels. The fully polarized 
stationary state is shown. 

Starting from equations (jHJ) we proceed in two steps which will be detailed in the following: 

1. We linearize these equations around the fully polarized steady state in which the only 
non-zero elements of p are p 44 = n and p o = N . 

2. From the linearized classical equations, interpreted as semiclassical equations for the 
mean vàlues of the collective operators, we derive the corresponding Heisenberg- 
Langevin equations. 
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B. Linearized Heisenberg-Langevin equations 



By linearization around the fully polarized solution we obtain equations for the "fluctua- 
tions" or deviations of the from their steady-state vàlues. Such linear equations coincide 
with the linearized semiclassical equations for the collective atòmic operators operators mean 
vàlues: 



Pki 
Pki 



(Sik) 
(Iik) 



kj = l 6 



k,l = 9,0 



(5) 
(6) 



where Ski = Y%=i \k)i(l\i for k, l — 1, . . . , 6 and Iti = HiLi \k)i(l\i for k, l — 9, are the collec- 
tive atòmic operators in the metastable and ground state, respectively. The corresponding 
linearized Heisenberg-Langevin equation for the operators are obtained by adding zero-mean 
valued fluctuating terms which are the Langevin forces. Denoting by f a the Langevin force 
for the operator a we get a closed set of equations: 



S21 
S32 

S43 
Io9 



— 7m<S'21 + /21 

2 / 7 



9 7m 



7/2 2 \ 

-g7m f 5*65 - yV3S 2 i - jS 32 ) + f 6 5 



Jr, 

1 

3 



— -S'43 H (S32 + 5*65) 



v/3 



7/A)9 + /43 



-37/J09 + 27 m - -S 65 + ^- S'43 + 



fi 



on 



(7) 
(8) 

(9) 
(10) 

(11) 



If a and (3 denote two system operators, (f a {t)fp{t!)) = D a p5(t — t') where D a ^ is the cor- 
responding coefficient of the diffusion matrix which can be calculated using the generalized 
Einstein relations llj] for an ensemble of uncorrelated àtoms. The non-zero coeficients are 



-043,34 — |7 m n, -Dq9,34 — -^43,90 



^l m n, D 



09,90 



2j m n. 



(12) 



Langevin forces are necessary to the consistency of the model. Otherwise, the non- 
Hamiltonian character of the exchange terms leads to a violation of the Heisenberg un- 
certainty relations. Physically, these forces originate from the fluctuating character of the 
ME col·lisions and their correlation time is the collision time, much shorter (~ 10~ 13 s) than 
all the times scales we are interested in. 
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C. Consequences of the Heisenberg-Langevin equations for ME col·lisions 



We notice that Eqs. ©-© for S21, S32, 5*65 form a closed subset of equations. This 
means that in the frequency domain each of these variables can be expressed as a linear 
combination of the Langevin forces /21, ^32, f&b- However, in the fully polarized limit we 
consider here, these Langevin forces do not contribute to the diffusion matrix. It follows 
that these variables do not contribute to the spin noise and can be neglected. One is then 
left with only two equations 

543 = -^S43 + -^7/-f09 + /43 (13) 

y/3 

l09 = -7/^09 + 7m — 5*43 + fo9 ■ (14) 

Let us introduce the transverse spin quadratures S x , S y 

S x = (S u + S i3 )/2, S y = i(S u -S 43 )/2 (15) 

(and similarly for the ground state spin transverse components I x , I y ) and let us assume 
that the ground state is initially squeezed, while the metastable àtoms are in a coherent spin 



state. Integrating O-jHI) with the initial conditions Aljj(0) = AZjJ(0)/(iV/4) = e~ 2r and 



AS 2 (0) = A5?(0)/(n/4) = 1 one finds the normalized steady state variances to be 



3nN 

' y " ^ " J (3n + A^) 2 

N 2 



AS v = 1 ~ [1 ~ e ' 2r ] I ,n 2 ( 16 ) 



^^-^-^WTW (17) 

Since n N (typically n/N ~ 10 -6 ), the ground state spin is still squeezed by approxi- 
mately the same factor e~ 2r , whereas the metastable àtoms squeezing is negligible (in n/N). 
By introducing the correlation functions Cs and Ci of two individual spins in the metastable 
and ground state respectively: 



AS 2 -1 A/2-1 , 

this simple calculation shows that ME col·lisions tend to equalize the correlation function 
(up to some numerical constant): Cs = 3Cj. If the ground state spin is squeezed, Cj has 
a negative value of order 1/N, corresponding to significant collective correlations for the 
iV-particle ensemble. However, as n N, this negative value of the correlation function in 



the metastable state is by far too small to induce squeezing into the n-particle metastable 
state, which would require Cs ~ —1/n. For e~ 2r = 1 we recover the coherent spin state with 
no correlation between the ground state and the metastable spins. 

Noise spectra can also be derived in a similar fashion. By defining the noise spectrum as 

S XiXj (u) = ídr e-^ix^x^r)) (19) 



where Xi, Xj are fluctuations of system operators and for the same initial conditions Ai^(O) = 

AJj(0)/(JV/4) = e~ 2r and ASj(O) = AS 2 (0)/(n/4) = 1 we get: 

n(Ne~ 2r + 3n)N 2 5(uj) 9 lexc nN 

Iy ' Iv[U) 2(N + 3n) 2 18^ 2 + 2(iV + 3n) 2 7 2 xc 1 j 

ç , s 3n(Ne~ 2r + 3n)n 2 5(uj) 3-f cxc nN 

S *> S * [UJ) 2(N + 3n) 2 18lu 2 + 2(N + 3n) 2 7 2 xc 1 ' 

The equal time correlations and (|17|) can be recovered from these formulas by integra- 

tion: 

(xiXj) = ^- J du S XiXj {uS) . (22) 

For an initial coherent spin state (e~ 2r = 1), the ME collision process does not change the 
collective spin variances, but it affects the noise spectra. The 5-shaped atòmic spectra of 
the two spins in absence of ME col·lisions acquire a width of order 7 ex c(^ + 3n), that is, of 
order j m . The contribution to the total variance of the "broad" part of the spectrum which 
is not sensitive to initial squeezing in the system, is large for the metastable state and small 
for the ground state. 



III. THE MODEL FOR SQUEEZING TRANSFER 

In figure |2] are represented the 3 He energy levels which are relevant for our squeezing 
transfer scheme. The àtoms interact with a coherent control field of Rabi frequency Q 
and frequency u\ that we treat classically, and a cavity field described by operators A and 
A^ . The field injected into the cavity, A in with frequency u 2 , is in an amplitude-squeezed 
vacuum state: (A in ) = and AX 2 n = e~ 2r , AY^ 2 = e 2r , where we have introduced the field 
quadratures 

X = A + A ] Y = i(A ] - A) . (23) 

The coherent field (7r-polarized) and the squeezed vacuum (<7~-polarized) are tuned to the 
blue side of the so-called C 9 transition (À = 1.08 /mi) from the F = 3/2 level of the 2 3 S' 
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metastable state to the 2 3 Pq state, the highest in energy of the 2 3 P multiplicity [15]. The 



atom-field Hamiltonian of the system is: 



H = H + h {fie-^Sra + S 82 ) + A(g A S 74 + g B S S3 ) + h.c) 



(24) 



where H describes the atom-field free evolution, gA,(B) = dA : (B)y2nu; 2 /hV are the coupling 
constants between the àtoms and the cavity field, V being the volume of the cavity mode and 



d>A,(B) the atòmic dipoles of the transitions 7 <-> 4, 



3). The system is initially prepared 



in the fully polarized state (ioo) = N and (S44) = n by preliminary optical pumping. 

Non-Hamiltonian contributions to the evolution of the system operators describe damping 
of the cavity mode, spontaneous emission from the excited state and the ME col·lisions 
described in detail in the previous section. Linearizing the equations in the rotating frame 




ME 
col·lisions 



FIG. 2: a) Metastable and excited sublevels of 3 He. Three coupling constants to the light are 
introduced. b) Squeezing transfer scheme using a control field íl, a squeezed vacuum field A and 
metastability exchange col·lisions. 9 and are the ground state sublevels. 

around the fully polarized state solution we obtain the following closed set of equations: 



S 2 \ = -{lm -i&i 2 )S 21 + ittS sl + /21 (25) 

Ssi = -(7 - i(Ai8 - 25 las ))S 81 + iüS 21 + fsi (26) 
2/7 \ 

5*32 = ~lm [~2 S32 + ^ S<21 + Se5 J + ^23^32 ~ Í^(Ss8 ~ S 72 ) + f 3 2 (27) 

S 72 = _( 7 _ Z (A 27 - 25 las ))S 72 - iíl(S 78 - S 32 ) + f 72 (28) 

5*43 = lm ^"^43 + — g-(532 + S 65 )j + ^7/I 9 + ^34^43 ~ ÍttS 47 + /tó (29) 

7/2 2 \ 

5*65 = ~g7m í Sq5 - ljVSS 21 - -S32) + ÍS^S 65 + / 65 (30) 



S 



where 



S47 

S38 
S78 
A)9 
À 



-(7 + iA 47 )S 47 - igAnA - ÍÜS43 + U7 
-(7 + i^x)Sas ~ iiï(S S2 ~ S 7S ) + fa» 



-(27 - i5 87 )S 78 - iQ(S 72 - S 38 ) + /' 



78 



-37//o9 + 27 m (s 32 - —565 + "^"^43 + "^~^ 21 



(k + iA c )A - ig B S 38 - Í9aS 4 7 + v2kA 



+ iS 90 I 09 + /, 



09 



(31) 
(32) 
(33) 
(34) 
(35) 



Ajj = (Ej — Ei) — lü2 

5ij = (Ej — Ei) — 5l as 
°las — ^1 — ^2 , 



(36) 
(37) 
(38) 



7 is the coherence decay rate due to spontaneous emission from the excited state and col·li- 
sions and we supposed Q to be real. The non-zero atòmic diffusion coeficients are 



-043,34 — —J m n, -D 4 3,90 — -^09,34 



2V3 



7 m n, -D 47i7 4 = 27n, D 



09,09 



2 7m n (39) 



We notice that metastable variables S21, S$i, S32, £72, Ses, «S38 and 5*78 form a closed subset 
of equations involving Langevin forces which do not give rise to non-zero diffusion coefficients 
in the fully polarized limit we consider here. Using the same argument as in section |Hj we 
deduce that these variables do not contribute to the spin noise and can be neglected. One 
is then left with only four relevant equations 



5*43 = 5*43 + — 7// 9 + ÍS34S43 - ÍÜS 47 + / 43 

S47 = -(7 + iA)S 47 - igAnA - iQS 43 + fw 

Í09 = g (-37/^09 + 7mV / 35 , 43 ) + 2<5// 9 + /09 

À = 



(k + iA c )A — igAS 47 + v2kA 



(40) 
(41) 
(42) 
(43) 



with A = A 47 and 5i = Sgo- 



IV. NUMERICAL RESULTS 



Equations (|ï*ü|) - ()4*3*|) can be used to find the variances of the metastable and ground 
state spin numerically. A typical result is displayed in figure El for which we assume that 
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a squeezed vacuum field with AX? n = 0.5 is injected into the cavity with the coherent 
control field in the squeezing-transfer configuration. In this figure AS 2 and AI 2 represent 

B [Tesla] 

10" 7 10" 6 10" 5 10" 4 10" 3 10" 2 

II I III 1 — I I 1 1 II 

1 .0 ~ ) CD CD CD - 

0.9- *S 
0.8- 

0.7-, 
0.6- ffi 

ES 

EB 

0.5-i ■ ■ — 

■ — i — i ■ — 

10" 3 10" 2 10" 1 10° 10 1 10 2 10 3 

r/v 

■m 

FIG. 3: Symbols: numerical calculations for spin variances in ground state (squares) and metastable 
state (circles), as a function of the ratio T/^f m (lower x-axis). The corresponding magnètic field 
needed to satisfy the resonance conditions (|46jl and (|47|) is shown in the upper rr-axis. Numerical 
vàlues of parameters are e~ 2r = 0.5, C = 500, k = IOO7, A = -2OOO7, 7 = 2 x 10 7 s~\ 7 m = 5 x 10 6 
s . The crossed squares correspond to a calculation including an extra relaxation rate 70 = 10 3 
s _1 for the metastable variables. The lines correspond to the analytical predictions (|ïïü|) and (|ïïTj) . 

the variances of S y and I y , both normalized to their coherent spin state vàlues. They are 
plotted as a function of the ratio T/j m , where T is the pumping parameter 

T = 7 3Ü 2 (l + (7)/A 2 , (44) 

and C = g 2 n/(l·í'-f) the cooperativity. It is precisely this ratio T/7 m which acts as a con- 
trol parameter to decide how the available squeezing of the field is shared between the 
metastable and the ground state spin. If T ^> 7 m , correlations are established among the 
metastable-state spins, the leakage of correlation towards the ground state being negligible. 
The metastable collective spin is squeezed while the ground state spin remains unsqueezed. 
In the opposite limit r 7 m , spin exchange is the dominant process for metastable àtoms; 
they transfer their correlations to the ground state which then becomes squeezed, while the 
metastable state remains unsqueezed. 
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In this plot we have chosen the best conditions for squeezing transfer: 

1. The metastable coherence S43 is resonantly excited by the two fields in a Raman 
configuration. By introducing the effective two-photon detuning for this coherence 

5 = 5 3A + n 2 /A (45) 

accounting for the light-shift of level 3, this condition reads 5 = 0, or 

(£ 4 - E 3 )/h + Q 2 /A = u 1 -u 2 . (46) 

2. The ground state coherence /09 should be resonantly excited by the metastable coher- 
ence (<5f = 0), i.e. 

(Eo - E 9 )/h = wi - w a . (47) 

In practice a magnètic field (shown as the upper x-axis) can be used to simultaneously 
fulfill ÍI46|) and (|47|) . When the resonance conditions are fulfilled the difference in the Larmor 
freqüències in the metastable and in the ground state is exactly compensated by the light- 
shift induced by the coherent control field. Choosing T = 0.l7 m as a working point, the 
required field is about B = 57 mG, corresponding to uj = 184 Hz. 

The vapor parameters in the figure correspond to a 1 torr sample at 300 K, with 7 m = 
5 x 10 6 s _1 and 7 = 2 x 10 7 s -1 , and a metastable atom density of 3.2 x 10 10 atoms/cm 3 which 
gives n/N = 10~ 6 . The symbols with a cross are a second calculation in which we added 
a finite relaxation rate in the metastable state 70, to account for the fact that metastable 
àtoms are destroyed as they reach the celi walls. We notice that only the ground state spin 
squeezing in the region r «C 7 m is affected. 

V. ANALYTICAL RESULTS 

In order to have a better physical insight it is possible to find simple analytical results 
within some reasonable approximation. By adiabatic elimination of the polarization S47 and 
the cavity field assuming 7,k> y m , 7/, one obtains 

S43 + (— + - - tó)S 43 = -^J 09 + / 4 a - -^fé7 + ^\J~ A in (48) 

^09 + (7/ - ÍSj)Ï09 = + /OB (49) 
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In deriving ()48j) we assumed a Raman configuration A > 7, y < 1 and that the cavity 
detuning exactly compensates the cavity field dephasing due to the àtoms: Ac = Ckj/A. 
From equation (jlHj) we see that (j m + T)/3 is the inverse of the characteristic time constant 
for the metastable coherence evolution. 



Resonant case 



If the resonance conditions (J4òj) and (J47|) are satisfied (5 = 5j = 0) and in the limit 

7/ <C r, 7 m , we can calculate the variances of the metastable and ground state spins 



N K 7m C_ _ _ 27 i 

n { r c 



^-T -fTtcTï' 1 -^ (50) 



^í-iV-—^ 1 -^ < 51 > 



which are plotted as full lines in figure 01 



B. Non-perfectly resonant case 

In order to test the robustness of our scheme, let us now concentrate on what happens if 
the resonance conditions (}4"ü)) and ()47|) are only approximatively satisfied. We will focus on 
the variance of the ground state spin coherence io9- 

By adiabatically eliminating the metastable coherence S43 one obtains 

j 09 + [T F + ib) 7 09 = /09 + lm f ~ í/43 - ^/ 47 + Í^JÏaJI (52) 

7 m + r - i35 \ A A ïk J 

The real part in the brackets 

r (7m + r) + ( 3» m 

(ira + T) 2 + (3Sf 

is the inverse of the effective time constant for the ground state coherence evolution which 
would also be the "writing" (or "reading") time of the quantum memory. T^ 1 = 2s in the 
example of figure El for T = 0.l7 m . It would be proportionally shortened by increasing the 
metastable àtoms density although Penning col·lisions prevent in practice metastable àtoms 
densities exceeding 10 10 -10 11 at/cm 2 . The imaginary part in the brackets 

b = -( lf ^ — =r + Si) (54) 
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is a light-shift "brought back" to the ground state, which is zero in the resonant case. 
Equation (J52J) can be used to calculate the best squeezing (optimized with respect to the 
transverse spin quadrature) of the ground state coherence: Al£ est = mm e A/f with I e = 



I x cos 9 + I v sin 9. We obtain 



N 



C 



where 



r + 7m + (35) 2 /r c + i 



1 — (e r + msinh(2r)) 



m = 1 — 



i + (P/T F y 



(55) 



(56) 
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FIG. 4: Normalized ground state spin variance corresponding quantum transfer efficiency r/j as 
a function of Si/Tf (while 5 = 0) (plots (a) and (b)), or 5/(F/3) (while 5i = 0) (plots (c) and 
(d). Symbols: numerical integration of equations í[25 [) -(|35 [) . Lines: analytical expression (|55|1. 
Parameters are the same as in figure|31and T = 0.l7 m . 

We show in figurelUthe effect of a frequency mismatch in on the normalized spin variance, 



Aly and the corresponding squeezing transfer efficiency r\i 



Vi 



1 - AP y 



1-e 



-2r 



(57) 
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In this example, a frequency mismatch of the order of T/3 in the metastable state or of 
the order of IV in the ground state affects the efficiency of the squeezing transfer. The 
condition for the ground state frequency matching (J47)) imposes stringent requirements on 
the homogeneity of the magnètic field. Because of the sinh(2r) in equation (J55j) . the larger 
the squeezing the worse are the consequences of a mismatch in the condition on 8j = on 
the ground state àtoms. Physically, if a significant dephasing between the squeezed field and 
the ground state coherence builds up during the squeezing transfer time, the àtoms will see 
an average between the squeezed and the anti-squeezed quadrature of the field noise. We 
can easily estimate the required magnètic field homogeneity as follows. Let us introduce the 
Larmor evolution freqüències in the metastable and ground states: in low field, hu a = fi a B 
(a=I,S) with fij/h = 3.24kHz/G and fig/h = 1.87MHz/G, and let AB be the maximum field 
difference with respect to the optimal value in the celi volume. For low field, the condition 
on AB to preserve the transfer efficiency reads [ijAB <C hTp. Since Íl 2 /A ~ — 
we get ï^:^3^^ < 1 or ; i n the regime V <C j m , 600^^ <C 1. With the parameters 
of figure El this gives a condition on the magnètic field inhomogeneity: AB / B < 4 x 10~ 4 . 
In figureElwe calculated the variance of the ground state spin as a function of T/7 m for an 
increasing inhomogeneity AB/B from zero (thick line) to 6x 1CT 4 . In practice a homogeneity 
of 100 ppm should be sufficient for the chosen parameters to guarantee that all àtoms will 
be squeezed. 

VI. OPTICAL READOUT 
A. Outgoing field squeezing 

n 

As briefly stressed in |3[ the squeezed fluctuations which are stored into the nuclear spins 
can be retrieved optically in the field exiting the cavity by using the reverse transfer process. 
Indeed, once the write sequence of the quantum memory has been completed, both the fields 
and the discharge can be switched off, leaving the àtoms in the fundamental state in a spin- 
squeezed state. After a variable storage time, switching back on the discharge and only 
the control field in the same configuration as for the writing phase (T <C 7 m ), will rapidly 
put a small fraction of àtoms in the metastable state and start the reverse transfer process 
from the fundamental àtoms to the field. The correlations in the ground state will slowly 
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FIG. 5: Normalized "best" variance of the ground state spin as a function of T/j m (lower x-axis) 
for an increasing inhomogeneity AB/B from zero to 6 x 10~ 4 by steps of 1 x 10 -4 . On the upper x- 
axis we show the corresponding homogeneous magnètic field needed to satisfy resonance conditions 
and (@7J). Numerical vàlues of parameters are e~ 2r = 0.5, C = 500, k = IOO7, A = -2OOO7, 
7 = 2 x 10 7 s -1 , 7 m = 5 x 10 6 s -1 , 70 = 0. 

transfer via the metastable state to the intracavity field. This will then result in squeezed 
fluctuations for the field exiting the cavity, which can be measured by homodyne detection. 

More quantitat i vely, if we still assume that the metastable spin observables and the 
intra-cavity field adiabatically follow the ground state spin observables and the evolution 
equations for the fluctuations of the squeezed component are in the resonant situation, we 
have 



Íy{t) = -T F I y {t)+{3X in {t)+jy{t) 

X out (t) = V2^X{t) - X in (t) 



(58) 
(59) 



with 



P 



7m + T 2A V « 



(60) 
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Denoting by t = the start of the readout sequence and by e~ 2r = Ai^(0)/(iV/4) the initial 
squeezing of the ground state nuclear spin, the two-time correlation function of the outgoing 
field amplitude quadrature can be obtained via (JoUj) after integration of 



C{t, t') = (X out (t)X out (t')) = 5(t - t') - 2Y F ·q I [l - e- 2r }e- TFÍ t +t,) (62) 

The 5-correlated term corresponds to the vacuum fluctuations contribution, whereas the 
second term corresponds to a transient squeezing for the outgoing field which is proportional 
to the initial atòmic squeezing. In (Jü2*j) . r\i designates the optimal quantum transfer efficiency 
in the ground state 

Ira C /ro\ 

7 m + 1 G + 1 

The ground state squeezing can be adequately measure d by homodyne detection using a 
temporally matched local oscillator as shown in Refs. [ÏIQ]. Using a local oscillator with 
envelope £(t) the normalized power measured by a Fourier-limited spectrum analyzer inte- 
grat ing over a time T is given by 

P(t) = — / p / dr dr'e-^-^£{r)£{r')C{r,r') (64) 

C[t) A J~ti/T 27T Jt Jt 

In order to measure the atòmic squeezing one has to maximize the temporal overlap between 
the local oscillator and the field radiated by the àtoms: £(t) oc e~ rpt . For such a local 
oscillator and for an integration time longer than the readout time T^ 1 the measured power 
can be written as the sum of a shot-noise term M and a time-dependent signal term S 
proportional to the initial squeezing: 

P(t) =M-S[\ - e ~ 2r }e- 2TFt (65) 

with S ~ TjjAÍ. The ground state nuclear spin fluctuations can therefore be measured 
optically with the same efficiency rjj as in the write sequence. However, because of the slow 
character of the correlation transfer process in the ground state the readout time is as long 
as the write time. As expected it is not possible to access the quantum memory faster during 
the readout than during the write phase. One could think of a faster readout method by 
transferring the fundamental àtoms fluctuations to the metastable àtoms and perform the 
optical readout in the regime T ^> 7 m . However, as we showed in section[Tïl starting with 
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a squeezed fundamental spin and first switching on the discharge (without the fields) will 
transfer very few correlations from the fundamental to the metastable àtoms and almost no 
squeezing will be retrieved in the field. 



VII. ENTANGLING TWO SEPARATE SAMPLES 

A direct and important extension of the previous results is that it is possible to transfer 
quantum correlations between different light beams to two spatially separated nuclear spins. 
If one disposes of EPR fields this allows to entangle two separate ensembles [l^ . Such EPR 
atòmic states are very useful for quantum information protocols involving the manipulation 
of continuous variable entanglement, such as atòmic teleportation for instance [Ïq| . 

Let us consider two identical ensembles 1 and 2 illuminated by EPR-correlated vacuum 
fields and coherent control fields fij (i = 1,2). Without loss of generality we assume 
symmetrical field correlations of the form 

A 2 X% = AVi? = cosh(2r) (i = 1, 2) (66) 
<4?4?> = -V$Y<Í>) = sinh(2r) (z ^ j), (67) 

i.e. that the amplitude quadratures are correlated and the phase quadratures anti-correlated: 
A 2 (x£ } - x£ 2) ) = A 2 (r£ } + Y^) = 2e- 2r . For perfect entanglement (r = oo) these EPR 
variances vanish. Both spins are initially prepared in a coherent spin state and we assume an 
equal incident power on both samples (ííj = f^)- Under the same adiabatic approximations 
as before, the fluctuations of the transverse spin components satisfy equation of the form 

L·i = —^fIxí — (Win + Ixi-, (68) 
Íy% = —^Flyi + I3X\ÏI + f yi , (69) 

(i = 1,2). Because of the linearity of the coupling in this regime, the EPR atòmic nuclear 
spin operators, I x \ + I x2 and I y i — I y2 , are clearly coupled to the EPR field operators 

J t {hl + h2) = -T F (Ixl + 1x2) ~ P(Y^ + Y^) + f xl + f x2 (70) 
f t (Iyl - Iy2) = -T F (I yl - I y2 ) + PÍX® ~ X^) + f yl - f y2 (71) 
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The amount of EPR-type correlations between the incident fields is given by the half-sum 
of the EPR variances 

1 



£ f = - 

s 2 



A 2 (X^ - X%) + A 2 (y« + Y®) = 2e-- (72) 



In the Gaussian approximation the entanglement between the nuclear spins can be evaluated 
using the same quantity (also normalized to 2) 

S T = ^ [A 2 (I xl + I x2 ) + A 2 (I yl - I y2 )] (73) 

It follows from (j70H71|) that the last two quantities are simply related by 

£ I = r ]l £ f + 2(l-r ]l ). (74) 

Like squeezing entanglement can also be in principle perfectly mapped onto the nuclear spins 
with an emciency r// (fH3j). close to unity in the regime r ^> -y m and C ^> 1. Let us introduce 
the correlation functions Cj 1 '^ of individual spins inside the ensemble i (i=l,2): 



C}W = ^_i (< =i )2 ) (75) 

and the correlation function Cj 1 '^ of two individual spins belonging to the different ensembles 
i and j: 

C f) = i^M [i^j = 1,2) (76) 

where the overline indicates the normalization of the correlation functions to N/A. In our 
case for r\i ~ 1 we get: 

Cr = C?» ^ (77) 

Cr-C?'^^l. (78) 

It is interesting to note that the two correlation functions C} 1 ' 1 ^ and Cj 1 ' 2 ^ become approxi- 
mately equal for a large entanglement e 2r 3> 1 so that an individual spin is about as much 
correlated with the other spins in its own ensemble as with the spins of the other ensemble. 



VIII. THE IMPERFECT POLARIZATION CASE 



The nuclear polarization of the sample is defined as 

p _ (joo) - (-^99) 
" </oo) + (J99) 
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(79) 



In practice polarization between 80% and 85% are currently achieved by optical pumping 
in dilute 3 He samples [so | . If the àtoms are prepared in a state which is not fully polarized 
- F / 1 - the situation is clearly more ccnplicated than we described in a and in the 
present paper. In particular, equations (|25)) - (j35j) and (J39j) obtained by linearization around 
the fully polarized state are no longer vàlid. We did not perform a complete analysis in the 
P 1 case. However, one can have a good idea of the result by using the simplified model 
of 3J which involves only two metastable sublevels (see figure|ïïj). As in section ITnl a Raman 
transition is driven by a coherent control field of Rabi frequency Q and a squeezed vacuum 
cavity field: 

H = H + h{nS 31 e- Wlt + gAS 32 + h.c.} . (80) 



In this toy-model the control field Q also acts as an optical pumping beam (able to transfer 
the àtoms from sublevel 9 to sublevel 0) and we introduce explicitly a relaxation in the 
ground state, so that P ^ 1 in steady state. Let us introduce for this model the rescaled 




ME 
col·lisions 



FIG. 6: Sublevels 1 and 2 are metastable, level 3 is the excited state, 9 and are the ground state 
sublevels. We include a relaxation process in the ground state so that P ^ 1 in steady state. 

coupling constant g, the atòmic one-photon detunings Ai and A, the two-photon detunings 
S and <5j, and two pumping parameters T p and V: 



A 1 = (Es-Ej/h-U!, 
5 = Ai- A + 



A 



A 2 



A = (E 3 - E 2 )/h - u 2 , 
Ó! = (E - E 9 )/h - (ut - uü 2 ) 

r = r p (i + c), 



(81) 
(82) 

(83) 
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where 7 is the optical coherence decay rate and C is the cooperativity parameter defined by 
equation (jHJ). For the atòmic operators we introduce S + = S21, S- = S12 



S x = (£_ + S+)/2, S y = i(S- - S+)/2, S z = (S 22 - S u )/2, (84) 

and similarly for the ground state operators. In the limit of large one photon detunings 
the excited state and the optical coherences can be adiabatically eliminated, yielding a set 
of equations similar to those of Ref. Q| with the addition of metastability exchange. By 
adiabatically eliminating the field (assumed to be resonant in the cavity) and for 5 = 0, 
Sj = 0, we obtain: 

S + = -{T p + 7m )5 + + 7//+ + 2igAS s + +fs + (85) 
S z = -(T p + lm )S z + lf I z + + i~g[A^S + - AS-} + f §z (86) 



A = -(K + iA c )A + igS + + y2nA in (87) 
/+ = -( 7/ + rO J + + lm S + + f I+ (88) 
Í z = -(7/ + Ti)/, + lm ~S z + f h (89) 

The semiclassical version of equations ([8'K |) -(f8ü ]) has a stationary solution (S + ) = = 
{A in ) = and with 

(S z ) = P = 7— p-^ j wTp S (90) 

7m 1 + ri(r p + 7m )/(r p7/ ) 

We will have in practice Ti <C 7/, meaning that the nuclear polarization in the metastable 
state P* = (S z )/(n/2) and the nuclear polarization in the ground state P = (I z )/(N/2) 
are almost equal. In this toy-model the stationary P is determined by the balance between 
the decay Ti and the pumping T p . In reality, the àtoms will be previously pumped more 
efficiently with resonant light. When we linearize the equations around the steady state we 
obtain 

= -(f + lm )S + + ^P*V2^A rn + 7/ J+ + f § (91) 
U = -(7/ + + lm S + + f I+ (92) 



with 



f = T p (l + CP*) (93) 
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Starting from equations (jïïT|) -(|ü2 ^) one can proceed as in section|V|to obtain 

AI 2 1 n' P* — 1 



where 



AI ^mr^p +{p ' e ~ lH+ é— (94) 



C = CP* and T f = ~ 7/F (95) 

r + 7 m 

J^- (96) 

c + i r + 7m r7 + r\ 



For P* ~ P and C ^> 1, we have finally 



A/2 = ^e- 2r + (l-^)/P (97) 

Equation (|97|) shows that the main consequence of having P ^ 1 is a rescaling of the 
cooperativity and the pumping parameter V and the quantum transfer efficiency rj'j, which 
are reduced by a factor P. Let us note that, for 1, when no squeezing enters the cavity, 
the àtoms are no longer in a coherent spin state. This shows, however, that strong squeezing 
transfer is still possible with a non-ideal polarization. 
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IX. APPENDIX 

Evolution equations of the density matrix elements under ME col·lisions are: 



d 




dt pu = 


Texc 


d 




dt pl2 = 


Texc 


d 




dt pl3 = 


Texc 


d 




dt pu = 


Texc 



N Pll + ^P99 (P22 + 3pn + 2p 55 )^ 

í -Np 12 + -p99 ({p23 + pm) V3 + 3pi 2 ) + -^P90 (P22 + 3pn + 2p 55 ) 
12 \ 

N Pl3 + gP99 (Pl3 + P24) + gP90 ((P23 + P5g) + 3p 12 ) ) 



\/3 N 

-Np U + ^~P90 (Pl3 + P24) 
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d 

dt p22 



dt p ™ 



dt pu 

d 

di p33 



d 

dt m 



d_ 

di 
d_ 

di 



dt 



d_ 
di 



Pm 



dt 
d_ 

di 
d_ 

di 



Poo 



P09 



P99 



7e 



7e 



Texc 
Texc 



Texc 

P44 = 7cxc 
P55 Texc 

P56 Texc 



Texc 
Texc 
Texc 



2 

-Np 22 + -P99 (2p 2 2 + P55 



P66 + 2p 3 3) + ^P90 (v / 3p21 + P65 + P32 

2 1 \ 

+ gP09 (-\/3pi2 + P23 + P5e) + gPoo (P22 + 3pn + 2p 55 ) J 

-iVp 2 3 + ^P99 (v / 3p34 + P56 + P23) + ^P90 (2p 2 2 + P55 + P66 + 2p 33 ) 

^ 2 / r \\ 

,Pl3 + P24) + -POO (V 3pi2 + P23 + Pm) I 

2 1 \ 

-ÍVP24 + -P90 (VS (p 23 + P56) + 3p 34 ) + gP00 (Pl3 + P24) j 

1 2 

-Np 33 + -P99 (2p 6 6 + P33 + 3p 44 ) + gP90 (v / 3p 4 3 + P65 + P32) 

2 2 \ 

+ -p 09 (v / 3p 3 4 + P56 + P23) + gPOO (2p 2 2 + p55 + Pm + 2p 3 3)j 
-iVp 3 4 + -g-P90 (2p 6 6 + P33 + 3p 44 ) + ~p00 (v^3 (p 23 + Pm) + 3p3ï) 



-Np AA + -poo (2p 66 + p 33 + 3pu)^j 

1 2 

-^p 55 + gP99 (2p 2 2 + P55 + Pm + 2p 3 3) - gP90 (v / 3p 2 l + P65 + P32) 

2 2 \ 

■-P09 (v^p^ + P23 + Pm) + gPoo (P22 + 3pn + 2p 55 )j 

2 1 

~Np 56 + -pgg (v / 3p 3 4 + P56 + P23) ~ gP90 (2p 2 2 + P55 + P66 + 2p 33 ) 

-^pP09 (Pl3 + P24) + ^Poo (VSp 12 + P23 + Pse)^ 
2 2 

-iVp 66 + -p 99 (2p 66 + p 33 + 3p 44 ) - -P90 (v / 3p 4 3 + P65 + P32 

2 1 

-gP09 (v / 3p 3 4 + P56 + P23) + gPOO (2p 2 2 + P55 + P66 + 2p 33 ) 
-"-Poo + g (3p 4 4 + P66 + P22 + 2p 55 + 2p 33 ) iV 
-™P09 + g# ((P43 + P2l) V 7 ^ + 2p 32 - p 65 )^) 
-WP99 + ^ (P33 + 2p 22 + 3pn + p 55 + 2p 66 ) iV 
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